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SUMMARY 


Studies have been conducted to determine the feasibility of using aileron 
or spoiler controls as alternates to pitch control for large horizontal 
axis wind turbines. The NASA Mod-0 100 kw machine was used as the basis 
for this study. Specific performance studies were conducted for 20% 
chord ailerons over the outboard 30% span, and for 10% chord spoilers 
over the same portion of the span. Both control systems utilized control 
deflections up to 60°. Results of the study show that either ailerons 
or spoilers can provide the control necessary to limit turbine power in 
high wind conditions. The aileron system, as designed, provides overspeed 
protection at hurricane wind speeds, low wind speed starting torque of 
778 N-m at 3.6 m/s, and a 1.3-1. 5% increase in annual energy compared to a 
fixed pitch rotor. The aileron control system preliminary design study 
includes aileron loads analysis and the design of a failsafe flyweight 
actuator for overspeed protection in the event of a hydraulic system 
fail ure. 
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INTRODUCTION 


Current large scale horizontal -axis wind turbines utilize variable pitch 
of all or a portion of blade span for starting, stopping, and modulation of 
speed. However, design of variable-pitch turbine blades of large size poses 
several structural and mechanical problems. Fixed-pitch rotors for large wind 
turbines offer advantages over variable pitch rotors or tip-controlled rotors. 
The blades of the fixed-pitch rotor are oriented so that the fatigue-producing 
cyclic gravity loads are in the direction of the greatest dimension (the sec- 
tion chordwise direction). This minimizing of cyclic stresses results in 
weight and cost savings. A fixed-pitch rotor does not require blade pitch 
bearings or the complex bearing assemblies and linkages needed for pitch con- 
trol and further stress reduction results from combining a teetering hub with 
fixed pitch rotor. These factors contribute to significant cost savings of 
fixed-pitch rotors over variable-pitch rotors. 

The fixed-pitch rotor, however, has two inherent problems: (1) start-up 
initiation and shutdown control, and (2) prevention of overspeed. Setting 
fixed-pitch blades at an angle to maximize annual energy may result in near- 
zero starting torques when the rotor axis is aligned with the wind, since the 
blades are completely stalled. It has been demonstrated that yaw control may 
be used for start-up of fixed-pitch wind turbines. However, yaw control is 
not fast enough to prevent overspeed. Yaw systems operate slowly-, as much as 
one minute may be required to yaw sufficiently to prevent overspeed. Typical 
overspeed rates are so high that the rotor speed can double in half that time. 

Although yaw control might be used to provide starting and to shut down 
at low wind speeds, other means are required for overspeed control. One pos- 
sible concept is the use of retractable spoilers on the outboard portion of 
the blades to reduce rotor speed by decreasing section lift and by increasing 
drag. An alternative concept is the use of hinged trai ling-edge surfaces 
(ailerons) which can act similarly to spoilers in preventing overspeed, and 
can also act to improve rotor torque at less-than-design wind speeds, thus 
improving annual energy output. 

Spoiler and aileron control systems have been studied under a 6-month 
grant from NASA Lewis Research Center. This report presents results of these 
studies which include aerodynamic performance of aileron and spoiler controls 
and their effects on turbine performance, selection of a control system sized 
for the NASA MOD-O 38.1 m diameter wind turbine, and preliminary design of 
mechanical systems for control surface actuation and fail-safe operation. 


PRELIMINARY CONSIDERATIONS 
Control System Options Studied 

Geometric variables available to the designer for the present study are: 
control surface type (aileron or spoiler), spanwise extent of control surface, 
control surface chord, and for the case of the spoiler, hingeline location. 
Most of the torque produced by a propeller- type wind turbine is produced by 
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the outer portion of the blade. For this reason, torque modulation devices 
should be located on this portion of the blade. The MOD-O turbine uses pitch 
control over the outer 35% of the blade. Control surface span for aileron and 
spoiler studies was selected as the outer 30% (see figure 1), 


Ailerons 


Ailerons have several characteristics which make them attractive as can- 
didates for wind turbine aerodynamic control. First, they produce nearly the 
same aerodynamic effects as spoilers of the same chord and deflection angle. 
Second, ailerons have the capability of providing positive as well as negative 
increments in lift. Furthermore, downward deflection of ailerons will increase 
maximum lift coefficient. The ability to increase lift at low wind speed con- 
ditions is attractive, since increased lift can be expected to result in in- 
creases in torque and hence power production. 

Spoiler Type 

Plate and plug type spoilers have both been successfully employed on air- 
craft (see sketch below). 




Plate-Type Spoiler 


Plug-Type Spoiler 


Sketch A - Spoiler Types 


The two types are essentially equal in terms of lift and drag increments pro- 
duced for a given projection height. The plug type have low actuation loads, 
but reduce torsional stiffness and strength of the airfoil more than the plate 
type. The plate type also causes less penalty if mis -rigged, and they have 
hinge moments which tend to close them. Based upon these considerations, plate 
type spoilers were selected for the present application. 

Spoiler Hingeline Location 

Experiments have shown that spoiler effectiveness increases as the loca- 
tion of the spoilers is moved forward (ref. 1). Unfortunately, forward spoiler 
locations result in a time lag between actuation and change in lift. For this 
reason, forward spoiler locations have rarely been selected for airplane con- 
trol systems. In the case of a wind turbine, time response requirements are 
not so severe. In aircraft, control response must occur within the time re- 
quired to travel a few chord lengths. Low solidity wind turbines such as the 
MOD-O machine have blade tips which travel a large number of chord lengths per 
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revolution (more tfian 100 chord lengths at mid-control position), and control 
response within one revolution is considered quite satisfactory. A second con- 
sideration for spoiler location is the effect on the structure of cutting into 
the airfoil, and a third consideration affecting the choice of spoiler type 
and location is the effect of "mis- trimmed" or "mis-rigged" spoilers on drag 
at conditions for which spoiler deflection should be zero. In addition, any 
gap or protuberance associated with the spoiler or actuating mechanism will add 
to the basic section drag, and the penalty associated with drag of this type is 
much more severe for forward locations than for aft positions. It is a cardinal 
rule of airfoil design that leading edges should be kept as smooth and unclut- 
tered as possible. Based upon these considerations, the present design studies 
were conducted for spoilers with hingeline locations aft of 50% chord. Because 
of the availability of data for spoilers with a 70% hingeline, that location 
was selected. 


PERFORMANCE ANALYSIS 


Leading-Edge Suction Studies 

Traditionally, aerodynamic characteristics of airfoil sections are stud- 
ied in terms of lift and drag coefficient components of the resultant force. 
Another set of components, sometimes used, are normal and chordwise force co- 
efficients. For a wind turbine blade element with zero twist and zero pitch, 
only the chordwise force produces torque. It is, therefore, instructive to 
study the airfoil characteristics in terms of this torque producing component, 
called "leading-edge suction." The relationships between c £ , cd* and the 
chordwise and normal coefficients are: 

c s = C£*sin(a) - Cd * cos (a) (1) 

c n = C£*cos(a) + c c j*sin(a) (2) 


where: 


c s = suction force coefficient, positive forward 
c n = normal force coefficient, positive upward 

Figure 2 shows the c^, Cj and c_ as functions of angle of attack for the 
NACA 23024 airfoil with 20% aileron. A theoretical relationship between c s 
and angle can be derived for a symmetrical section from thin airfoil theory, 
as follows: 


c^ = 2irsin(a) , Cd = 0 (3) 

Substituting these value relationships into equation (1) above, the following 
result is obtained: 

C s = 2irsin 2 (a) (4) 

In real flow, drag is present, and at zero angle of attack c s is equal to -Cd- 
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At angles of attack near stall, the suction decreases abruptly, and again 
becomes negative. Thus the useful angle of attack range for developing positive 
power is from the angle at which enough "theoretical" suction is developed to 
overcome the basic section drag, up to the angle of attack at which stalling 
occurs. For the NACA 23024 airfoil the range of angles for positive suction is 
between 2° and 17°, A second region of positive suction appears at about 60° 
and remains up to 90°. 

The relationship between wind speed, rpm, and blade angle is shown in 
figure 3. When the turbine axis is aligned with the wind, the angle of attack 
on the blade will range from 90° at zero rpm to 0° at infinite rpm. For an 
untwisted blade element at zero pitch angle, the angle of attack will be equal 
to the angle <i>. For low lift coefficients and low blade solidity, the induced 
velocity factors a and a' will be small. It is useful to assume the limiting 
case of zero induced effects to study trends. Setting a= a' = 0, 

tan(a) = tan(<j>) = w ~ (5) 

0)1 


Single Point Performance Analysis 

Selecting the 75% radius station as a representative station for blade 
element analysis, this gives: 


tana = 


Vi nd 
oj( . 75 R) 


( 6 ) 


Rearrange: 


^wind 

.75 R (tan (a)) 


(7) 


Using this relationship, and a specified wind speed, the rpm can be calculated 
corresponding to a given angle of attack. From the c s curve for the NASA 23024 
airfoil, it is seen that zero torque (zero c s ) will occur at several angles. 

For a wind speed of 54 m/s (120 mph), equilibrium rpm values are calculated as 
shown in Table 1: 


Table 1 - 

Equilibrium States at Zero Torque 

Angle of Attack Wind Speed 

RPM 

2° 

54 m/s (120 mph) 

1029 

17° 

54 m/s (120 mph) 

116 

60° 

54 m/s (120 mph) 

21 


This study shows that the highest rotational speed corresponds to the lowest 
angle of attack. 
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The question of which of the possible equilibrium states will really exist 
Is answered by considerations of stability. Examining the graph of c s vs. a, 
and realizing that torque Increases asr c s Increases, and that rpm decreases as 
a increases, the question of stability Is easily understood. 

Sketch B shows these characteristics: 



Sketch B - Stability 


Since an increase in rpm from point 1 will result in negative torque, and 
a decrease in rpm will result in positive torque, it is clear that this point 
is in stable equilibrium. Point 2, on the other hand, has the opposite charac- 
teristics. An increase in rpm from point 2 will result in an increase in 
torque, accelerating the rotor toward point 1, and a decrease in rpm from 
point 2 will result in negative torque, decelerating the rotor toward point 3. 
Thus, point 2 is unstable . These same considerations will show that point 3 
is a second stable equilibrium point. 

These arguments are easily generalized to state that equilibrium points 
with positive slope on the graph illustrated will be stable, while equilibrium 
points with negative slope will be unstable. To protect the wind turbine from 
overspeeding in high wind with loss of electrical load on the generator, it is 
necessary that either no stable equilibrium states exist, or that the stable 
equilibrium "runaway" rpm be low enough that no damage will be sustained, and 
the remaining rotor energy can be absorbed by the system mechanical brake. 
Sketch B shows that point 1 is the most dangerous equilibrium point, since it 
corresponds to the highest stable equilibrium rpm. 

Overspeed Control 

In the event of electrical failure of a generator being powered by a wind 
turbine, it is necessary to provide means for preventing overspeed or "runaway" 
which would result in overstressing of the rotor structure and subsequent fail- 
ure. For large wind turbines this poses a difficult problem for a mechanical 
brake design. The preferred method is to use aerodynamic braking such as pitch 
control or similar means. The braking device must be capable of functioning 
not only at the rated or design wind speed, but also at the highest wind likely 
to be encountered. For the present study a hurricane wind speed of 54 m/s 
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020 mph) was selected as the maximum wind speed. Combining the requirements 
of a "double emergency" situation consisting of generator failure along with 
hurricane wind results in a design requirement to prevent overspeed with zero 
torque and wind velocity of 54 m/s . 

Figure 2 shows the relationship of c s and angle of attack for aileron de- 
flection of -60°, and spoiler deflection of 60®. The results for 60° upward de- 
flection of either spoiler or aileron show that the region of positive cs in the 

low (pre-stall) angle of attack range has been eliminated. Thus, no positive 
torque will be developed over this range of angles. At angles of attack beyond 
60° however, positive cs persists even with aileron or spoiler control surface 
deflected upward 60°, and positive torque will be developed. From the argu- 
ments given earlier, an equilibrium "runaway" state can be expected at condi- 
tions corresponding to 60° angle of attack. For the hurricane wind speed of 
54 m/s, the equilibrium rotational speed will be about 21 rpm. Since this ro- 
tational speed is well below the operating speed of 33 rpm, it does not pose 
a hazard in terms of centrifugal blade loads. Mechanical braking must be used 
to bring the rotor to rest. 

From these preliminary studies based upon simplified analysis and hand 
calculations, it is indicated that either a spoiler or an aileron system will 
be adequate for providing overspeed protection to the turbine. The next sec- 
tion documents more sophisticated computer studies of turbine performance with 
aileron and spoiler controls, which confirm these preliminary findings. 


Computer Analysis of Performance 

Analysis of performance of the MOD-O wind turbine was conducted using the 
PROP computer code (see ref. 2), modified as explained in Appendix A. The pro- 
gram receives input information about rotor geometry, control mode (tip section 
pitch, spoilers, or ailerons), and operation (varying rpm or wind speed) through 
a series of data cards. For the analysis reported here, the following input 
items were considered constant and were not varied: 


B = 2 
BO = .97 
XETA= 0.167 
HB = 2 ft 
R = 62. 5 ft 
THETP = 0 deg 
THETI = Odeg 
H= 100 ft 
SI = 0 deg 


(number of blades) 

(NASA tip loss model) 

(velocity power law exponent) 

(hub radius) 

(turbine radius) 

(pitch angle) 

(twist angle at all inboard stations) 
(hub altitude) 

(coning angle) 


Blade planform geometry used in all calculations corresponds to the wooden 
inboard blade (shown in figure 1). The length of the spoilers and ailerons is 
30% of radius. This geometry was entered by means of data cards. 
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Input data which were varied during this study included: 

V freestream wind velocity 

OMEGA angular velocity of turbine rotor 

MCON method of control (tip section pitch, spoilers, 
or ailerons) 

TIPICH angle of pitch of tip section 

AILDEF angle of deflection of ailerons 

SPDEFL angle of deflection of spoilers 

The program permits use of various sizes of radial increments in the 
calculation of performance. Early in the analysis studies, comparison was 
made of results of computation using radial steps ranging from 10% to 2%. To 
minimize error, 2% increments were used in all subsequent calculations. 

In operation, the computer program starts at the outboard end of the 
blades and calculates components of forces for each incremental radial blade 
element, working inboard to the hub. 

Local chords and amount of twist are determined by the subroutine SEARCH 
using linear interpolation. Subroutine CALC is then called to calculate 
axial and rotational interference factors (a and a 1 ), local angle of attack 
and lift and drag coefficients. Relative wind relationship to the blade sec- 
tion is shown in figure 3. The relationships between angle of attack and 
lift and drag coefficients of the NACA 23024 airfoil are obtained from sub- 
routine NACAXX and modified by subroutines TIPLOS (tip and hub losses) and 
INCREM (increments of c^ and cj due to aileron or spoiler deployment). 

As can be seen in figure 3, angle of attack is a function of a, and a is 
a function of the section lift. Thus, evaluation of a is an iterative process 
which is discussed in detail in Appendix A. The program may be controlled to 
limit a between -0.5 and +0.5, or to set a = 0, or to permit a to be unlimited. 
Following computer runs to study the effectiveness of these schemes, remaining 
performance predictions were calculated with the restriction that -0.5<a<0.5. 
Radial interference (a 1 ) was assumed zero for all the present studies. Local 
elements of force in the plane of rotation (producing torque) and normal to 
the rotor plane (producing thrust) are determined by iteration at both ends 
and at the center of each blade section, and the forces are integrated across 
the blade to give total quantities for the blades. Details of the airfoil 
aileron, and spoiler modeling are given in Appendix A. 

Power Modulation at High Mind Speeds 

For wind speeds greater than the wind speed for which the electrical gener- 
ator is rated, it is necessary to either take the wind turbine off-line or to 
"dump" power by aerodynamic or other means. With an aileron or spoiler control 
system this "dumping" is accomplished by deflecting the control upward to re- 
duce lift and torque. Figure 4 shows power coefficient versus tip speed ratio 
for the basic rotor, and the rotor with full up aileron, and with full up 
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spoiler. These results also show that either the aileron or spoiler can pro- 
vide the control necessary to reduce power to zero over the normal range of tip 
speed ratios. Because of the potential of the aileron to provide starting 
torque and to provide increased power at low wind speeds, the aileron was se- 
lected for further studies. 

The data in figure 4 show that the aileron chord could be reduced consi- 
derably, perhaps as much as 50%, and still provide necessary overspeed protec- 
tion. While smaller aileron chord would reduce cost, it was decided to retain 
the 20% aileron chord size for two reasons. First, if one aileron became inop- 
erative the rotor could be stopped with the remaining aileron. Second, it is 
known from experimental tests that the PROP code and similar theoretical methods 
under- predict maximum power by a considerable amount. Therefore, a substantial 
margin must be provided in the design of a control system using these analyti- 
cal methods. If experimental evaluation of the new control system reveals that 
it is oversized, corrective action can be taken at that time. Improved aero- 
dynamic modeling for wind turbines should be developed to reduce discrepancies 
between theory and experiment in order to reduce costs of experiments and 
avoid expensive design errors. 

Figure 5 shows how power modulation would be accomplished with the 20% 
aileron control system for a synchronous generator operating at 33 rpm. For 
wind speeds lower than 8.2 m/s (18 mph) with zero control deflection, the 
power developed is lower than the rated value of 120 kW shaft power, corres- 
ponding to 100 kW of electrical power. For wind speeds in excess of 8.2 m/s 
(18 mph), the turbine is capable of producing more power than the generator 
can absorb. By deflecting the control surface upward, power can be modulated 
to provide the required 120 kW at speeds up to 12.9 m/s (29 mph). For speeds 
greater than this value, it will be possible to provide 120 kW by using down 
aileron deflection, an option not possible with spoilers. Thus the aileron 
system can provide rated power up to 15.6 m/s (35 mph). 

Overspeed Control 

Computer studies with full up aileron were conducted at the hurricane 
wind speed (54 m/s) condition described earlier, for rotational speeds ranging 
from 10 to 30 rpm. These runs show that the zero torque "runaway condition" 
occurs at 19 rpm. This value compares very favorably with estimate of 21 rpm 
from the simplified analysis given in Table 1, and confirms that the proposed 
control system will provide adequate overspeed protection. 

The control system must be capable of providing actuation forces suffi- 
cient to overcome the closing aerodynamic hinge moment. The critical design 
condition is maximum overspeed rotational speed of 40 rpm with 60° up aileron 
at zero wind speed. The hinge moment for this case is 186 N-m per aileron. 


Aerodynamic Starting 

Wind turbines with pitch control are ordinarily started by rotating the 
blades to an unstalled angle, and then decreasing the pitch angle as rpm 
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increases. With a fixed pitch rotor, very little starting torque is developed, 
and other means for starting must be employed. 

Recent tests by NASA with the MOD-O turbine have demonstrated that start- 
ing can be effected in a fixed pitch mode by yawing the rotor away from the 
wind. While this mode of operation is possible for a large scale machine, it 
is inconvenient in that yaw rates must be low to prevent over-stressing the 
rotor system. 

With an aileron control system it is possible to generate starting torque 
without yawing the rotor away from the wind. At zero rpm the effective angle 
of attack is 90°, which means that the basic blade will be fully stalled, and 
torque will be nearly zero. With up aileron, however, a larger starting torque 
is provided. The computer code used for performance studies is not capable of 
calculation at zero rpm, so runs were made at very low rpm values, and the re- 
sults were extrapolated to zero rpm. This study shows that the optimum aileron 
angle for starting is -60°, and that a starting torque of 778 N-m (574 ft-lb) 
is available at a wind speed of 3.6 m/s (8 mph). Because of the dynamic pres- 
sure effect, this torque will be proportional to the square of the wind speed. 
Table 2 shows starting torque for several wind speeds. 

Whether starting by use of ailerons can actually be effected depends upon 
two factors wich have not yet been determined. First, the friction of the 
generator and gear train must be overcome. Second, once starting is achieved, 
rpm must be increased to the point that the blades become unstalled so that the 
operating speed of 33 rpm can be achieved. Accelerating the rotor to an un- 
stalled state cannot be simulated with the present computer codes. Full-scale 
rotor tests will be required to evaluate the practicality of this starting 
method. 


Table 2 - Starting Torque with 

60° Up Aileron 

Wind Speed 

Torque 

2.2 m/s ( 5 mph) 

304 N-m 

3.6 m/s ( 8 mph) 

778 N-m 

4.5 m/s (10 mph) 

1216 N-m 


Power Increase at Low Wind Speeds 

Another interesting control option is possible with the aileron system. 

At wind speeds lower than rated, positive aileron deflection can be utilized to 
increase power available. 

In order to assess the benefits of aileron deflection at speeds below 
rated, it is necessary to calculate annual energy output for competitive control 
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schemes using some mathematical model for the wind probability distribution. 
In the present study, a Rayleigh distribution has been used. Table 3 shows 
results of a study of this type, comparing zero control with optimum aileron 
deflections. Optimum aileron deflections are 5° or lower. 


Table 3 - Use of Ailerons 

to Increase Annual 

Energy 

Mean Wind Speed 

5.36 m/s 

5.81 m/s 

6.26 m/s 

(12 mph) 

(13 mph) 

(14 mph) 

Increase in , R0/ 

Annual Energy "° /o 

1.6% 

1.3% 


These results show that using optimum aileron provides small gains in 
annual energy. If the rated power were higher, or mean wind speeds lower, 
larger increases would result. 


PRELIMINARY DESIGN 


General 

The general configuration of the aileron control system (shown in figure 
6) was selected based upon airplane flight control design experience, two- 
dimensional wind tunnel research of aileron and spoiler control systems, and 
performance studies of the MOD-O wind turbine using the PROP computer program, 
with the airfoil subroutine modified to accommodate aileron control character- 
istics. From these latter studies, it was determined that an aileron of 20% 
chord and 60° upward deflection having a span extending over the outer 30% of 
the rotor would provide adequate overspeed protection. The 15° down travel is 
adequate for providing added power at low wind speed conditions. A bell crank 
transfers spanwise motion of a primary push rod to chordwise motion. The pri- 
mary push rod is actuated by a hydraulic cylinder located in the inboard por- 

tion of the rotor. 

The counter-weight shown in the figure acts under the influence of centri- 
fugal force to drive the aileron to the full up position, providing fail-safe 
operation in the event of hydraulic system failure. The hydraulic cylinder 
works against the weight to bring the aileron to a desired setting. The deci- 
sion to locate the hydraulic cylinder on the rotating blade was made largely 
because of existing hydraulic hardware used for pitch control on the MOD-O 
machine. In other applications it would be quite feasible to locate the hy- 
draulic cylinder within the rotor hub, and transfer the needed spanwise actua- 
tor motion through a push rod extending to the rotor hub. A mechanical latch 

locks the ailerons in the full up position until released by restoration of 

hydraulic system pressure. 
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It is desirable to locate the aileron hingeline near the airfoil upper 
surface rather than on the mean line for two reasons. First, when deflections 
of as much as 60° are to be used, the airfoil surface forward of the aileron 
must be cut back to a position well forward of the hingeline to clear the de- 
flected aileron. This clearance produces a large gap in the airfoil surface 
when the aileron is in neutral position; this gap is an added source of drag. 
Moving the hinge toward the upper surface minimizes the gap size and the drag. 
Second, upper surface hinging permits maximum moment arm inside the airfoil con 
tour for attachment of the aileron actuating mechanism. 

On a straight-tapered blade of constant percentage thickness most para- 
meters that need to be calculated can be stated in terms of local blade chord 
which, in turn, depends upon blade station (radial distance from the center of 
rotation). 

The Dimensional quantities given in this section of the report are in 
English units rather than S.I. units because of current fabrication practice. 

In the calculations which follow, two aileron planform arrangements are consi- 
dered. For one planform the aileron extends from the 69.1% radius (518.34 in. 
station) to the rotor tip, and for the second case the aileron extends from 
64.7% radius (485.00 in. station) to the rotor tip. Performance studies show 
that either arrangement will be satisfactory. The present design loads analy- 
sis will also indicate that no critical loading situation will occur for either 
geometric arrangement. For the present planform geometry, local blade chord is 
given by: 

c = 7.05786 - .0867925 x (8) 


where c = local chord, ft 

x = local radial distance, ft 

Airload Hinge Moment 

Data from several sources indicate that the maximum value of closing hinge 
moment coefficient, c^, that can be expected with the aileron deflected up (down 
wind) 60° will be less than 0.58. To be conservative, a value of 0.58 will be 
used for the present design study. 

Maximum hinge moment will occur when the angle of attack of the blade is 
nearly 0°, i.e., when rotational velocity is high and wind velocity is low. In 
this case, resultant velocity on the blade outboard sections will come almost 
entirely from rotational velocity and can be taken as v = xw. A rotational velo 
city of 40 rpm (4.189 radians/sec) will be used for design. 

Measured along the upper surface from the hingeline to the trailing edge, 
the aileron chord is .1930 of the local chord c. Thus, at any blade station, 
the incremental airload hinge moment will be 


H 


/ c h 


q c 2 dx 


(9) 
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where H = hinge moment 
q ■ (wx) 2 

c a = .193c 

and c is given by equation (8). 

Using sea-level standard air density and 40 rpm, the following results 
are obtained: 

H Inboard Aileron Station 

2412 in-lb 518.34 in 

2808 in- lb 485.00 in 


Aileron Weight 

The flat-pattern area of sheet metal required to make the aileron skin at 
any station is approximately .5557c per unit span. Total skin weight is cal- 
culated by integrating across the span of the aileron: 


r dA s 

w s = J w*t*-j^-*dx 


( 10 ) 


where w = specific weight of aluminum, 0.10 lb/in 3 

t = skin thickness, .020 in 
dA s . 

= skin area per unit span = .557c 

c = a linear function of x, from equation (8) 

The sheet metal area required for a rib at any station is approximately 
.765c for flanges plus .01240c 2 for the rib web. Treating ribs as distributed 
weight, rib weight equals 


f dA r 

w r = J w*t*-^-*dx 
dA r 

where = rib area per unit span - rib area/rib spacing 
dA r , 

-jjpr = (.765c + .0120c 2 )/rib spacing 

(The rib spacing in this case is 10 inches.) 


( 11 ) 
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Hinge and trailing edge strip will contribute about .027 Ib/in. Rivets, 
etc., are estimated to weigh about 1.0 lb total. The total weight of one ail- 
eron is given by: 


w total " w skin + w ribs + W hinge + TE + w rivets 


( 12 ) 


Total Weight Inboard Aileron Station 

17.8 lb 518.34 in 

21,01b 485.00 in 

Total Centrifugal Force on Aileron 

The total centrifugal force on the aileron is obtained by a span-wise 
integration: 


F c = /“ 2x ^ dx O 3 ) 

Hm 

where ^ = mass per unit span - weight per unit span/g 

Using the weight equations developed earlier for various components of the 
aileron, and assuming "distributed" rib mass as earlier, the following result 
is obtained for w=40 rpm. 


Centrifugal Force Inboard Aileron Station 

511 lb 518.34 in 

582 lb 485.00 in 


Centrifugal Moment about Aileron Hinge Line 

Since centrifugal forces on incremental mass elements of the aileron act 
neither through, nor parallel to, the aileron hinge line, these forces produce 
moments about the hinge line as shown in Sketch C on the following page. 

The hinge moment at any station is the local centrifugal force times the 
sine of the angle between the hinge line and a true radius, multiplied by a 
moment arm equal to the perpendicular distance between the hinge line and the 
line of action of the centrifugal force through the local section c.g. The 
magnitude of these moments are calculated below. Hinge mass is neglected since 
the moment arm is so short. 

For ribs, skin and trailing edge strip, the c.g. can be taken as being 
.053c below (upwind of) the hinge line. This is the dimension "d" in Sketch C, 
The direction of the hinge line, relative to the center-line of the steel tube 
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(a) Rotor and Centrifugal Force 


Sketch C - Centrifugal Moment 

spar of the blade is .05066 radians. The direction of the centrifugal force at 
the c.g. of any aileron section relative to the spar center-line is approximate- 
ly .6152 (c/r) radians. Thus e= .05066 + .6152(c/r) radians. 

Multiplying the force component by its moment arm about the hinge line 
and integrating over the aileron span gives: 

H c - /(» 2 )(x)(sin E )(d)|>dx (14) 

e is a small angle s so e*sine. 
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Substituting all values and integrating, the following results are ob- 
tained: 


Centrifugal Hinge Moment Inboard Aileron Station 

68.4 in-lb 518.34 in 

75.2 in-lb 485.00 in 

These values are much lower than the aerodynamic hinge moment. 

Design of Centrifugal Force Drive System 

Maximum airload hinge moment to be overcome = 2808 in-lb . To assure ex- 
cess capacity to drive the aileron to full up position, design for 2808x 1.5 = 
4212 in-lb moment. The centrifugal moment produced by the aileron mass will be 
neglected for conservatism in this analysis. 

Push rod moment arm about hinge line at full up aileron = 3.05 in. Thus 
push rod load = 4212/3.05 = 1321 lbs. Push rod moment arm about bell crank 
pivot is also 3.05 in. 

Counter-balance mass c.g. has a moment arm of approximately 11 inches when 
aileron is 60° up. Thus, the centrifugal force needed to supply 4212 in-lb 
moment is 4212/11 = 383 lbs . 

The mass is located at x=510 in. Therefore the centrifugal load factor 
at 40 rpm (4.189 rad/ sec) will be: 


Substituting, 



05 ) 


n = 23.18 


the weight required will be 


w = 


F 

n 


w = 16.52 lb 


Thus a centrifugal weight of this size will be able to overcome the aero- 
dynamic hinge moment at 40 rpm and 60° up aileron with a 50% margin of safety, 
and taking no credit for the centrifugal effect of the aileron mass. 

Sizing Hydraulic Cylinders 

Moment arm of the actuating cylinder force to oppose the centrifugal force 
will be approximately 3.9 inches with aileron up 60°. Assuming (in the extreme) 
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that a) is still 40 rpm and airload hinge moment has fallen to zero, the actuat- 
ing cylinder load required to start the aileron down will be (4212/3.9) = 

1080 lbs . 

If it is desired to produce this force with hydraulic pressure of 1000 psi 
or less, net effective piston area of the activating cylinder needs to be: 

A = £ (16) 

A = 1.08 in 2 


Assuming a piston rod diameter of 0.50 in, a piston diameter of 1.375 in 
will be more than adequate. 

Total piston travel available should be 5.5 in. The fluid volume change 
associated with full travel will be: 


Vol = J (D q 2 - D-j 2 ) (L) 
Vol = 7.1 in 3 


(17) 


Coordinating Aileron Deflection on the Two Blades 

Because of physical constraints imposed by the configuration of existing 
MOD-O rotor design, it is desirable that position coordination between the 
ailerons on the two blades be accomplished hydraulically rather than by mech- 
anical interconnect between the two ailerons as originally proposed. Assuming 
identical actuating cylinders and aileron drive linkages in the two blades, 
this requires some means of guaranteeing equal fluid volume flows to or from 
the two actuating cylinders as aileron position is changed. 

Since it is unlikely that loads in the two linkages will always be the 
same even though aileron positions may be the same, maintaining equal pressures 
in the two actuating cylinders will not satisfy the equal volume change require- 
ment. Thus, merely splitting the flow from one pressure source will not do. 

One simple way to guarantee equal volumes is to use some such arrangement 
as that shown on Sketch D on the following page. This arrangement provides 
a closed hydraulic system for each blade. Two identical master cylinders 
(shown as two chambers in a single housing) are driven by a single primary 
actuator. The two master cylinders and the cylinder that drives them consti- 
tute the unit labeled "motion equalizer assembly" on the drawing. 

Except for the effects of leakage, such an arrangement will keep positions 
of the two slave cylinders coordinated even though pressures in the two systems 
are different. Since no hydraulic system can be absolutely leak-free, it seems 
advisable to provide a means of supplying makeup fluid to each of the two 
closed systems to compensate for small amounts of leakage. Such a system would 
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require position sensors on each of the slave cylinders and on the motion 
equalizer assembly. When either slave cylinder position does not agree with 
equalizer unit position within acceptable limits, an injector (perhaps solenoid 
actuator) would inject a small amount of fluid into that system. It is empha- 
sized that the proposed arrangement would compensate only for very low leakage 
rates. Any gross leak would result in shutdown for repairs. 


Equalizing Twin Master Actuating Cylinders 



Blade # 1 
aileron 


Blade # 2 
aileron 


All of the items described above could be mounted on the rotor. If this 
were done and the primary hydraulic supply were mounted on non-rotating struc- 
ture, one fluid transfer gland on the rotor shaft and slip rings to supply power 
for the injector systems would be required in addition to whatever slip rings 
were needed for instrumentation and control circuits. 

The primary supply system could also be rotor-mounted thus eliminating all 
fluid transfer glands. However, the number of slip rings probably would be con- 
siderably increased if this were done. It appears that means could be provided 
to mount both the motion equalizer package and the pressure supply package on 
the outside surface of the rotor hub. At 40 rpm and a radius of 2,0 ft, centri- 
fugal loads will be only 1.1 g. It should be possible to mount all required 
hardware well within a 2.0 ft radius. 

Irrespective of where the components are located, the hydraulic cylinders 
will be sized as follows. The full area of the pistons in the twin master cylin- 
der will be available. Thus the diameter will be: 



where A = 1.08 in 2 . 

Thus 

D = 1.17 in 


Pressure in the single cylinder driving the twin master cylinders will be 
essentially the same as in each of the twd driven cylinders if the single cylin- 
der is made 1.75 inches in diameter. Its total volume flow requirement would 
be 5.5 x 1 .75 2 ir/4 = 13.23 in 3 . 
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To allow the ailerons to move quickly from normal operating position to 
full up-position for emergency shutdown of the wind turbine s the line(s) from 
the dump valve to reservoir must be designed to allow a high flow rate at low 
pressure. 

Aileron Twist 

Since the aileron actuation will take place at the inboard station, and 
the aerodynamic hinge moment is distributed along the span, the aileron will 
be subjected to torsional loading. The deflection in torsion (twist) is cal- 
culated by analyzing the aileron as a thin-walled cylinder (non-circular). The 
appropriate equation (from ref. 7) is: 

6 ■ Iwr^ (,7 > 

where T = applied torque per unit span 

z = perimeter of thin-walled member (aileron) 

A = cross-sectional area of aileron 
G = elastic modulus in shear of aluminum (3.8 xlO 6 psi) 
t = skin thickness (.020 in) 

T, z, and A are all functions of local chord, which in turn depends on x. This 
analysis ignores rib and hinge contributions to torsional stiffness. The design 
condition is again the case of 40 rpmwith 60° up aileron. The results of this 
analysis are given below: 

Total Twist Aileron Inboard Station 

3.30° 518.34 in 

3.88° 485.00 in 

This amount of twist is considered to be quite acceptable. 

Aileron Latching 

A spring-loaded latching system is proposed to hold the ailerons in full 
up position once they have been fully deflected. Hydraulic pressure would 
work against the spring for the subsequent re-start. This could be accomplished 
by providing for over-travel where the actuator cylinder is attached to the cen- 
trifugal mass arm so that initial piston motion when hydraulic pressure is re- 
applied will release the latch and further motion will bring the aileron toward 
the neutral position. 
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On sudden loss of electrical load, a dump valve in the pressure supply 
system would be opened. Overspeed and/or manual shutdown signals also would 
be fed into the control for this valve. The same signals would be used 
to interrupt power to the pump drive motor. 

Design Summary 

The design information is summarized as follows: 


Aileron Material: ,020 in thick 2024T4 aluminum 

(ribs and skin) 

Aileron Inboard Station 



518.34 in 

458.00 in 

Aileron Weight: 
(each) 

17.8 lb 

21.0 lb 

Centrifugal Force: 
(radial ) 

511 lb 

582 lb 

Centrifugal Hinge Moment: 
(opening) 

68 in-lb 

74 in-lb 

Aerodynamic Hinge Moment: 
(closing) 

2412 in-lb 

2808 in-lb 

Aileron Twist Under Maximum Moment: 

3.3° 

3.9° 

Centrifugal Weight: 

16.5 lb 

19.2 lb 

Hydraulics (1000 psi system) 



Actuating Cylinder 
Actuator Force: 

1080 lb 

1257 lb 

Piston Diameter: 
(w/0.5 in. D. shaft) 

1.375 in* 

1.375 in 

Travel (minimum): 
Flow Equalizer 

5.5 in 

5.5 in 

Twin Piston Diameters: 

1.17 in 

1.17 in 

Driving Piston Diameter: 

1.75 in 

1.75 in 


This cylinder is slightly oversized to be a standard dimension. 
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CONCLUSIONS 


1. Either ailerons or spoilers could be used to prevent runaway and 
to modulate power at higher than rated wind speeds. 

2. An aileron system was selected as prime candidate for the present 
applications based upon considerations of providing power output gains at 
low wind speeds, and the possibility of starting without yawing the rotor 
away from the wind. 

3. Preliminary loads analysis including hinge moments have been 
developed for a 20% chord, 30% span aileron control system for the MOD-O 
turbine. 

4. The system, as designed, provides overspeed protection at hurri- 
cane wind speeds, low wind speed starting torque of 778 N-m at 3.6 m/s, 
and 1.3 to 1.5% increase in annual energy compared to a fixed-pitch rotor. 


RECOMMENDATIONS 


1. Aileron hardware should be developed for testing the aileron control 
concept on the MOD-O turbine. 

2. Theoretical models for wind turbine performance should be improved 
to diminish the disparity between theoretical and experimental power values. 

3. Wind tunnel tests should be conducted to provide a better data 
base for airfoil characteristics at high (up to 90°) angle of attack, with 
and without control surface deployed. 
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Figure 1 - Blade Planform Geometry 
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Figure 3 - Blade Element Velocities 
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Figure 5 - Effects of Aileron Control on Power at 33 RPM. 
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APPENDIX A 


MODIFIED COMPUTER CODE: PROP 


Basic Program, "PROP" 

The computer program used in the studies conducted under this NASA grant 
is PROP . This program was originally written by Stel Walker and Robert 
Wilson at Oregon State University for NASA (ref. 2). Prior to being supplied 
to Wichita State University, some modifications were made to the program at 
Lewis Research Center to facilitate its use on the NASA computing system. 

After the program was received at WSU, it was modified to permit opera- 
tion on the university IBM system 370 computer. The program consists of a 
main program and six subroutines: 

Main Program: PROP 

Subroutines: SEARCH 

CALCZ 
TIPLOS 
NACAXX 
INCREM 
BESSEL 

As explained in the Performance Analysis section of this report, forces 
are determined at local stations and then summed for the complete blade. It 
was determined that it was desirable to use 21 radial increments. The compu- 
ter program starts at the outboard end of the blades and calculates components 
of forces for each incremental radial blade element, working inboard to the 
hub. 


Local chord and amount of twist are determined by the subroutine SEARCH 
using linear interpolation. Subroutine CALC is then called to calculate the 
axial and rotational interference factors (a and a 1 ), the local angle of attack 
and lift and drag coefficients. The relationships between angle of attack and 
lift and drag coefficients are obtained from subroutine NACAXX and modified 
by subroutines TIPLOS (tip and hub losses) and INCREM (increments of c % and 
Cd due to aileron or spoiler deployment). 

Convergence Problems and Related Modifications 


As can be seen in figure 3, angle of attack is a function of a, and a is 
a function of the section lift. Thus evaluation of a is an iterative process 
involving the following steps: 

1. Assume initial values of a = 0, a' = 0. 


2. Calculate <j> = tan" 1 ^ + (Al) 
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3. Calculate a = * - 6 


(A2) 


4. Calculate c^ and c^ as functions of a, including any increments 
due to spoiler or aileron deflection. 


5. Calculate c y = c £ cos<|> + c<jSin<f> (A3) 

and c x = c^sin* - c^cos* (A4) 

c y 

6. Calculate a = nf= (A5) 

8 sin 2 * c y 

— c x 

7. Calculate a' = (A6) 

8 sin * cos <j> 


8. Compare a or a' to previous values to check for convergence. 
If not converged, return to step 2 and repeat. 


This process had been programmed in subroutine CALC to perform the itera- 
tion up to 40 times with damping applied on the fourth, tenth, and fifteenth 
cycles. It was found that at tip speed ratios, X, of more than 15, the aver- 
age values of a were predicted to be greater than 0.5 for the straight blade 
at zero pitch angle. For blades with negative or positive loading outboard— 
as with tip pitched, ailerons deflected, or spoilers deployed— high values of 
a were obtained at lower values of X, with local values of a=1.0 (or greater) 
being obtained. While these unrealistic values were usually obtained at low 
wind speeds outside the normal operating range of the rotors, it was felt that 
they cast doubt on all performance predictions. 

In order to improve calculation of a, the possible number of iterations 
was increased to 50, damping was applied to each cycle, and local values of a 
were limited to -0.5<a<0.5. Also, an optional procedure was made available 
which disregarded the induced velocity by setting a = 0. Figure A-l compares 
the results of calculations with (1) a unlimited, (2) -0.5<a<0.5, and (3) a = 0. 
For the cases shown here, methods (1) and (2) give the same Cp. 

Elements of force in the plane of rotation (producing torque) and normal 
to the rotor plane (producing thrust) are determined by iteration at both ends 
and at the center of 2 % span blade section. The forces are integrated across 
the section using Simpson's rule and then simply summed to give total quanti- 
ties for the blades. 
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Airfoil Characteristics Subroutine, "NACAXX" 


Airfoil aerodynamic characteristics used in the PROP computer code are 
made available to the computer in the form of subroutine NACAXX. This tech- 
nique permits the program user to develop the aerodynamic information in a 
form most suited to his own needs, or most compatible with available data. 

The task of the subroutine is to determine lift and drag coefficients for a 
given angle of attack, and to return these coefficients to the main program. 

In subroutine NACAXX used for the present performance studies of the NASA MOD-O 
wind turbine, the aerodynamic characteristics of the NACA 23024 airfoil are 
stored as a set of algebraic equations, with logic to branch to the equation 
appropriate for a particular angle of attack range. The angle of attack at a 
given station is calculated from wind speed, angular speed, radius, and in- 
duced factor, a. 

In order to analyze an airfoil with a control surface, it is necessary 
to identify geometric characteristics of the control surface, such as control 
surface chord and deflection, as well as airfoil angle of attack. In addi- 
tion, it is desirable to store information concerning control surface hinge 
moments so that control surface actuation loads may be calculated. In the pre- 
sent research, it was decided that control surface data should be treated as 
increments to be applied to the basic airfoil characteristics. There were two 
reasons for adopting this procedure: first, no data were available for the 

NACA 23024 section with spoilers and ailerons; and second, a control surface 
incremental subroutine would have application to a variety of airfoils with 
minimal change. 

Incremental Subroutine, "INCREM" 


This subroutine receives control surface type, control surface chord, 
and deflection angle as input information. Output information is incremental 
lift coefficient, incremental drag coefficient, and control surface hinge mo- 
ment coefficient. The incremental data are added to the basic airfoil coeffi- 
cients, so that net or total coefficients are returned to the main program for 
evaluating blade element performance. 

The incremental data used in this subroutine were obtained from wind tun- 
nel experiments conducted on a 21% thick airfoil at WSU, specifically the 
LS ( 1 ) -0421 airfoil as documented in reference 3. Figures A-2 and A-3 show 
incremental data from the reference for lift and drag coefficients, and con- 
trol surface hinge moment data as a function of control surface deflection and 
angle of attack. These incremental effects should be reasonable for the NACA 
23024 and other airfoils of similar thickness. For sections less than 15% 
thick, more appropriate data should be used. 

The forms selected for the equations used to fit the incremental control 
surface data have been designed to have rather broad applicability, in terms 
of control surface chord, even though the present data are obtained from 
experiments with 20% chord ailerons and 10% chord spoiler, respectively. The 
resulting equations are believed to be applicable to control surfaces with 
chord lengths differing by as much as 50% from the source data values, but 
could lead to serious errors for larger chord differences. 
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Furthermore, the spoiler data are from experiments with a spoiler hinge- 
line location of 70% chord. While using the present data base to approximate 
spoiler characteristics for hingelines from 60% to 90% chord would probably 
not result in large errors, hingeline locations forward of 50% would be ex- 
pected to produce large errors, since it is known that spoilers are much more 
effective at forward locations. 

Within these constraints, the forms chosen for the incremental equations 
will have reasonable applicability. The coefficients should be adjusted based 
on more appropriate experimental data, when such data are available. 


Aileron Increments 


Lift: 


c a 


A = —(11.5 sin 6 a - 5 sin 2 6 a ) 

0° < <5 a < 60° , ! a l - “stall 

Ac* = — • (11 .5 sin 6 a + 3.35 sin 2 5 g ) ] 
-60° < fi a < 0°, i“l- “stall 


(A7) 


(A3) 


Drag : 


Ac d = ^ (1.05 sin 2 s a ) 

0°<6 a <60°, { ct | < “ s t a n 


Ac d * (°- 65 sin2 5 a) 

-60° < 6 a < 0° , I ct | < a s tan 


(A9) 


(A10) 


“stall = 
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Hinqe Moment: 


e h = ~ (-.0035 - ,0542a - 2.00 sin 6 a ) 
0°<6 a <60°, |a| <o $tall 

c h “ (-.0035 - .0542a - 2.95 sin 6 a ) 

-60° < 6 a < 0° , |a|<a stall 


(All) 


(A12) 


To determine overall hinge moment required to actuate a control surfaces 
integration of local section hinge moments across the span of the control is 
accomplished in program subroutine CALC, in parallel with other spanwise in- 
tegrations of blade torque and thrust loads. 

Post-Stall Behavior of Airfoil with Control Surface 


Post-stall characteristics of airfoils have been investigated for only 
a few section shapes. This lack of data is unfortunate, since fixed pitch 
wind turbines operate over an angle of attack range from approximately 0° 
to 90°. The aerodynamic characteristics of the NACA 23024 airfoil as modeled 
in subroutine NACAXX for angles of attack beyond stall were based upon wind 
tunnel tests of the symmetric 12% thick NACA 0012 section. For post-stall 
angles of attack with spoiler, it is assumed that the spoiler effectiveness 
vanishes and that the basic section data apply. This assumption is believed 
to be reasonable since for high angles of attack, the airfoil will separate 
in front of the spoiler, rendering it nearly ineffective. 

For an airfoil with aileron deflected, on the other hand, the lower sur- 
face shape is changed and the lift and drag will be different, even for 90° 
angle of attack. Since no experimental data for airfoils with ailerons at 
extreme angles of attack are available, a model was developed for the present 
study based upon theoretical results from potential flow analysis of flat 
plates using free streamline techniques. Theoretical results for Kirchoff 
flow, modified by correcting to a more realistic base pressure are used to 
develop aileron loads for post-stall cases (see ref. A2j see Sketch A2, p. 37). 

The normal force for this flow is given by: 

N = c n y P V 2 c (A13) 
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where 


_ 27rsin(a) 
n 4 + it s i n (a ) 


(A14) 


For a = 90° , equation (A14) gives Cp = 0.88. From two-dimensional experi- 
ments c n «2. The discrepancy between theory is primarily because the base 
pressure is not modeled properly by the Kirchoff theory. To correct for this 
effect, a correction factor K is introduced. 


c 


n 


2irsin(a) 

4 + ir | sin(g) 


(A15) 


where 

K = M=2. 27 


This equation (A15) will be used to calculate the normal force for a de- 
flected aileron at angles of attack beyond stall, where no experimental data 
are available. The load distribution will be assumed uniform over the aileron 
The loading on the portion of the airfoil forward of the aileron will be cal- 
culated in the usual way, using data from SUBROUTINE NACAXX. The steps in 
the total calculation process are: 

1. Obtain cz and c^ as a function of g for the airfoil without 
aileron from NACAXX. 

2. Transform lift and drag coefficients to chordwise and normal 
force coefficients: 
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= C d * cos(o) - * sin(a) 

(A16) 

= * COS (a) + C d * sin(a) 

(A17) 


3. Calculate the contributions of the area forward of the aileron: 


C C 1 3 c c* 0 - c a/c) 

(A18) 

c ni 3 c n * (1 - c a /c) 

(A 19) 

Calculate the aileron normal force coefficient 
into components parallel and perpendicular to 

, and resolve 
the airfoil chord: 

r „ 2.27 * 2it * sin(a+ 6 a ) 

L n2 — 

4 + ir * (sin (a + 6 a ) | 

(A 20) 

normal component: 


Ac n a = cn 2 *^-* cos(6 A ) 

(A21) 

chordwise component: 


ACc a = c n 2 *T* S1n ( 6 a^ 

(A 22 ) 


5. Calculate aileron hinge moment: 

c n? 

c h — 

(Assumes uniform load distribution on aileron.) 


(A23) 
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6. Calculate airfoil plus aileron force coefficients: 




c n = c n] + acn a 

(A24) 



C C 3 C C] + ACc a 

(A 2$ 

7. 

Convert 

to lift and drag coefficients: 




c l 3 c n *cos(o) - c c *sin(a) 

(A26) 



c d = c c * cos (a) + c n *sin(a) 

(A 27) 


{These 

are the net section lift and drag coefficients, including 


aileron 

effects . ) 


8. 

Return 

to main program. 




Spoiler Increments 


Lift: 


Ac* a (-9.9 sin 6 S ) 




0° < 6 S < 60°, I a | < <* sta1] 

' (A28) 



AC * 3 o, H>cc stall 

(A 2^ 

Drag: 


c<- " 

Ac rf = — (0.93 sin 2 6_ ) 
u e s 

(A 30) 

► 



°° - 6 s - 60 °» l°l - “stall j 
AC d = 0, N>a stal1 

(A 31) 
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Hinge Moment: 



(0,175a - 6.7) (sin 2 <5 S ) 
°° ioJ “stall 


! 


c } 

c h = (0.075a - 6.7) (sin 2 6 $ ) 

‘“stall so<0 ° 
c h ’ °’ l“l >a stall 


(A 32) 


(A 3 3) 


(A 34) 
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Expanded Output Modifications 


In order to facilitate analysis of effects of spoiler or aileron deflec- 
tions on wind turbine performance, the output of the computer program was 
expanded. 

The computer program output was modified to provide two alternative types 
of display of the performance results. The "long print-out" shows local val- 
ues of velocity, angle of pitch, angle of attack, coefficients of lift and 
drag, and the resulting local contributions to torque, thrust, and a. Inte- 
grated results are also printed including torque, thrust, power, c p , average 
a, and hinge moment. The "short print-out" displays only the integrated 
values. 

The reason for the long print-out is that it permits an investigator to 
examine local wind angles and angles of attack across the radius from hub to 
tip, the resulting blade loading and induced effects. 

An additional optional print-out was added to the program during the pre- 
sent research. At each radial station, the computer seeks to satisfy contin- 
uity and momentum relationships by an iterative procedure involving a, the 
induced velocity factor, net angle of attack on the blade section, and the 
resulting section lift and drag coefficients. In order to examine this itera- 
tive process, an optional command causes the printer to list each successive 
estimate of local a, c^, c<j, and the factor a. If a does not converge in the 
50 iterations a warning is printed on the output sheet. 

Hinge Moment 

Hinge moments for the spoiler studies were obtained directly from the 
wind tunnel data of reference A3. For the aileron studies, however, hinge 
moment data from this reference are not directly applicable, because of the 
large difference in airfoil camber between the two sections, especially near 
the trailing edge. The equation for hinge moment coefficient is given by: 

c h = c ho + c h a a + c h 6 6 (A35) 


where: 

c^ * hinge moment coefficient 

c^q - hinge moment coefficient at a=0° and 6 = 0° 

(a function of camber) 

Ch a * rate of change of hinge moment with alpha, obtained from 
the data of reference 5 

c^ = rate of change of hinge moment with aileron deflection, 

0 from reference 5 
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Since the value of cho depends on camber of the section, and no experi- 
mental data were available for the NACA 23024 section, a value for this para- 
meter was obtained by integration of the pressure distribution over the air- 
foil trailing edge region using theoretical data from reference A4. 

Computer Program Listing 

A listing of the computer program with all the modifications described 
above follows: 
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C PROGRAM PROP MODIFIED TO SON PERFORMANCE STUDIES BY VARYING 

C MIND SPEED OR ROTATIONAL SPEED., 

C 

C 

G 

C THIS PROGRAM HAS DEVELOPED BY STEL N. HALKER AND 

C ROBERT E, WILSON,- DEPARTMENT OF MECHANICAL ENGINEERING, 

C OREGON STATE UNIVERSITY- ALL QUESTIONS OR PROBLEMS 

C CONCERNING THE PROGRAM SHOULD BE DIRECTED TO THE ABOVE. 

C 

C MODIFIED BY MEL SNYDER, BILL HENTZ, AND CYRUS OSTOMARI OF THE HIND 

C ENERGY LAB AND THE AIRFOIL RESEARCH GROUP AT WICHITA STATE UNIVERS 

C 

C - - - MODIFICATION DATES — JUNE TO OCTOBER 1979 - 

C 

C ....... MAIN PROGRAM 

C PROP CALCULATES THE THEORETICAL PERFORMANCE PARAMETERS OF A 

C PROPELLER TYPE HIND TURBINE, IT UTILIZES A SIMPSON'S-RULE 

C METHOD / THREE PASS TECHNIQUE OF NUMERICAL INTEGRATION. 

C 

C METHOD OF CONTROL MAY BE - - 

C (1) PITCHING OF THE TIP SECTION 

C (2) USE OF OUTBOARD SPOILERS 

C (3) USE OF AILERONS 

C 
C 

1 DIMENSION CHD (2,50) ,THETX (2,50) ,CIM(50) 

COMMON /HOPE/ R, DR,HB, B, V, X,THETP, AM OD, H, SI, GO, OMEGA, RHO , VIS, HL, 

1 PI,RX,W , NPBOF, APF, RADEG,COSSI,CS2,NF,RR (50) ,CI (50) ,THETI (50) , 

2T1, T2,T3,T4,T5,T6,T7,T8,XETA,BH,BO,CAI (50) ,CSI(50) , DAI (50) ,DSI (50) 
3 , MCON ,HMF, KFLAG 
DIMENSION PHD (2) ,PFACT{2) 

DATA KW , PHD, PFACT/1 , * HP KW % 1.341,1./ 

........ DEFINE INPUT AND OUTPUT INDICES...... 


KR=5 

KP=6 

IOUT=6 


...... ..READ INPUT DATA., ...... . 

READ (KR, 42) MODE, NINC, INCVEL, INCRPM, I PRINT, KFLAG 

READ (KR, 10) R,DR,HB,THETP 

READ (KR, 10) B,V, OMEGA, AMOD 

READ (KR, 30) H , SI , NF,GO, NPROF 

READ (KR, 40) HL, BO, APF,XETA, HH 

HEAD (KR, 4 1 ) -AFl , AF2 , AF3, AF4 , AF5, MCON, TIPICH # TIPL, SPAN, DEFL, PCTSPN 
READ (KR, 20) (£R(I) ,CI (I) ,THETI(I) ,CAI (I) ,CSI (I) ,-DAI (I) ,DSI (I) ,1=1 , 
*NF) 

WRITE (KP,60) 

WRITE (KP,61) 

RM=0. 3048*R 

WRITE (KP, 62) AFl, AF2, AF3,AF4,AF5, B,RM , SI 
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WHITE (KP,63) 

GO TO (84*85,86) ,HCON 

84 WRITE (KP, 64) TIPICH,TIPL 
GO TO 89 

85 WRITE (KP»65) SPAN, PCTSPN,DEFL 
WRITE (KP ,67) 

GO TO 89 

86 WRITE (KP, 66) SPAN, PCTSPN,DEFL 
WRITE (KP, 68) 

89 GO TO (19,21,22) ,KFLAG 
19 WRITE (KP, 78) 

GO TO 91 

21 WRITE (KP, 79) 

GO TO 91 

22 WRITE (KP,83) 

91 IF (IPRINT.EQ, 1) GO TO 90 
WRITE (KP, 77) 

WRITE (KP ,70) 

WRITE (KP,74) 

WRITE (KP , 75) 

WRITE (KP,76) 

WRITE (KP ,70) 

90 DO 500 NO=1,NINC 

IF (IPRINT.EQ. 2) GO TO 13 
WRITE (KP, 69) 

WRITE (KP,7 1) 

WRITE (KP ,7 2) 

WRITE (KP,73) 

WRITE (KP, 69) 

13 IF(MODE.EQ. 1)GO TO 6 
IF (MODE-2) 499,7,499 

6 VDEL=INCVEL 
V=V+VDEL 

GO TO 8 

7 RPMI NC=INCRPM 
OMEGA=OMEGA+RPMINC 

8 PI=3„ 1 41 593 
RADEG = PI/180. 

FPSMPH = 88 „/60» 

ISTOP=0 

SCALE=1. 

CONTINUE 

R = R*SCALE 
HB = HB* SCALE 
OMEGA = OMEGA/SCALE 
DO 4 1=1, NF 

CIM (I)=0«3048*SCALE*CI (I) 

4 Cl (I) = Cl (I) *SCALE 

IREAD=0 

STORE INITIAL VALUES 

51 = R 

52 = DR 

53 = HB 

54 = THETP 
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55 = V 

56 = OMEGA 

57 = SI 

C SI IS CONOING ANGLE (PSI) 

COSSI=COS (SI*RAD£G) 

cs2=cossi*eossi 

W=R*COSSI*OMEGA* PI/30. 

x=w/v/fpsmph 

INITIALIZATION AND CONSTANT PARAMETER CALCULATIONS ...... 

THETP=-THETP*RADEG 
V=V*FPSMPH 
OMEGA=OMEGA*PI/30. 

SI=SI*RADEG 
SUM 1 = 0.0 
SUM2 = 0.0 
QX=0 . 0 
QY=0 . 0 
TX=0.0 
TY=0 . 0 
PX=0.0 
PY=0 . 0 
AK= 1 . 

ASTOP=0. 0 

INDUCED VELOCITY FACTORS ARE SET EQUAL TO ZERO TO BEGIN ITERATION 
A=0. 0 
AP=0 . 0 
HMOM=0.0 
CONTRL=3. 0 

•CONTROL* COUNTS STEP IN 3-STEP SIMPSON RULE INTEGRATION 
RHO=0. 0023769 199*EXP(-0.297*H/1 0000. ) 

RX=R 

CORRECT FOR CONING 
R=R*COSSI 
REF =BO*R 
DR=DR*R 
DRO=DR 
HB=HB*COSSI 

' NN • IS NUMBER OF STEPS ACROSS SPAN, DETERMINED FROM INPUT VALUE 
OF DR 

NN= (R-HB) /DR+3. 

IF (IREAD. GT. 0) GO TO 5 
DO 11 1=1, NF 
RR(I)=R*RR(I)/100. 

THETI (I) =THETI (I) *RADEG 
CONTINUE 

SET R LOCAL EQUAL TO RMAX 
RL=R 


NUMERICAL INTEGRATION FROM TIP TO HUB ...... 

CAT=1. 

IF (GO.EQ. 2. ) CAT=2. 
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CLF A= 1 . 

IF (GO.LT. 2. ) GO TO 563 

IF (GO . EQ. 3, ) CLF A=0. 0 

IF (IREAD.GT.O) GO TO 211 

CALL SEARCH (RL,CX1 ,THETX1 , CA,CS, DA, DS) 

211 C=CX1 

THET = THETX1 

CALL CALC (RL, C, THET, QX ,TX, A , AP , F, CLF A , CAT, CL, CD, ALFA, CH , CA ,CS ,DA , 
2 DS) 

C RE- INITIALIZE *A* 

563 A=0.0 

CAT=0. 0 

C BEGIN MARCHING ACROSS SPAN FROM TIP TO HUB 

DO 100 L= 1 , NN 

C DO NOT MARCH INSIDE HUB 

IF ( ( RL-H B) , GE. DR) GO TO 50 
ASTOP=ASTOP+1. 

IF(ASTOP. GE.2.) GO TO 93 
DR= (RL-HB) 

50 IF (GO.LT. 3. ) GO TO 311 

C TEST FOR STEP IN 3-POINT SIMPSON 

IF(CONTRL.EQ.O.O) GO TO 311 
TIP=RL-DR 

IF (TIP. GT. REF) GO TO 312 
IF (CONTRL. EQ. 2. ) GO TO 311 
DR= (RL-REF) 

CLFO= (REF-TIP) / (RL-TIP) 

CLF=. 5*CLFO 
CONTRL=1- 
GO TO 311 
312 CLF=0. 0 

311 DR 2= DR/2 «, 

DT6=DR/(6.*COSSI) 

C INCREMENT 'B' BY »DR/2* 

RL=RL-DR2 

IF (CONTRL. EQ. 0.0) CLF=1. 

IF (CONTRL. EQ. 2.0) CLF= (CLFO+1 . ) /2 . 

IF (IREAD.GT.O) GO TO 202 
C FIND LOCAL CHORDS AND ANGLES 

CALL SEARCH (RL» CHD ( 1, L) ,THET X ( 1 ,L) , CA, CS,DA ,DS) 

202 C=CHD (1 ,L) 

THET = THETX (1,L) 

C CALCULATE VALDES AT R+DR/2 STATION 

CALL CALC (HL,C,THET,QXP1,TXP1,A, AP,F, CLF, CAT,CL,CD, ALFA,CH,CA,CS, 
2DA, DS) 

C INCREMENT R BY DR/2 AGAIN (TOTAL INCREMENT = DR) 

RL=RL~DR2 

IF (CONTRL. EQ. 0.0) CLF=1. 

IF (CONTRL. EQ. 1.0) CLF=CLFO 
IF (CONTRL. EQ. 2.0) CLF=1.0 
IF (IREAD.GT.O) GO TO 203 
C FIND LOCAL CHORDS AND ANGLES 

CALL SEARCH (RL, CHD (2, L) , THETX (2 ,L) , CA, CS, DA , DS) 

203 C=CHD (2,L) 
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THET = THETX(2,L) 

C CALCULATE 7 ALOES AT R*DR STATION 

CALL CALC {RL,C,THET,QXP,TXP,A,AP,F,CLF,CAT,CL,CD, ALFA, CH, CA, CS„ DA, 
2DS) 

C THETA=LOCAL BLADE TWIST+PITCH SETTING 

THETA=THET+THETP 
PCRL= { 1 00. ) *RL/R 
RLM=0.3048*RL 
CM=0.3048*C 

C NEXT 4 STEPS INTEGRATE TORQUE, THRUST AND POSER, USING 3-POINT 

C SIMPSONS RULE 

C 3 VALUES ARE AT R, R+DR/2,AND B+DR 

QYX=DT6* <QX+4.*QXPHQXP) 

QY=QY + QYX 

TY=TY+DT6* (TX+4. *TXP1+TXP) 

PY=PY+OMEGA*QYX 
HMOH=HMOM+HMF*DR/COSSI 
IF (CONTRL. EQ. 2. ) CONTRL=0.0 
IF (CONTRL.EQ. 0.0) GO TO 313 
IF{ (RL-TIP) .EQ.0.0) GO TO 313 
IF {CONTRL. EQ. 1.) DR=REF-TIP 
IF (CONTRL. EQ. 1.) CONTRL=2. 

GO TO 314 

313 DR=DRO 

314 CONTINUE 
QX=QXP 
TX=TXP 

C NEXT 3 STEPS CALCULATE TERMS TO EVALUATE AVERAGE A (* AVA') 

C OVER THE SPAN 

C ‘AVA 'STEP BELOW COMPLETES THE CALCULATION 

AREA = DR* (2, *RL+DR) 

SUM 1 = SUM 1+AREA 
SUM2 = SUH2+ ( 1 .—A) * AREA 
IF (IPRINT. EQ. 2) GO TO 100 
29 GO TO (33,31,31) ,MCON 
33 CH=0. 

31 WMS=W*0. 447040 

WRITE (KP, 80) PCRL,RLM, CM,THETA, ALFA,W MS,CL,CD,CH,QXP,A, AP 
100 CONTINUE 
93 CONTINUE 

CTY=TY/(.5*RHO*V**2*PI*RX**2) 

CPY=PY/ (. 5*RHO*V**3*PI*RX**2) 

TP=PY/737.6 
PWR=TP*PFACT (KW) 

AVA= 1. -SUM2/SUM 1 

C RESTORE INITIAL VALUES FOR POSSIBLE RERUN 
IF(IREAD. EQ. 1) GO TO 105 
DO 104 1=1, NF 
RR (I) =1 00. *RR (I) /R 

104 THETI (I) =THETI (D/RADEG 
IREAD=1 

105 CONTINUE 
R = SI 
DR = S2 
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HB = S3 
THE TP = S4 
? = S5 
OMEGA = S6 
SI = S7 
S8 = X 

VMS=0.447040*V 
QYNH=QY*1» 355833 
TYN=4„ 44827*TY 
CQ=CPY/X 
HMOM=HMOM*» 7376 
IP (IPRINT. EQ. 2) GO TO 490 
WRITE ?KP,77) 

WRITE (KP, 7 0) 

WRITE ?KP,74) 

WRITE (KP, 75) 

WRITE ?KP,76) 

WRITE ?KP, 7 0) 

49 0 WRITE ?KP,81) V, VMS, OMEGA, X, THETP, QYM , TP, CPY,TYN,CQ, HMOM, AVA 
IF(IPRINT. EQ. 2) GO TO 500 
WRITE (KP,813) 

500 CONTINUE 
GO TO 502 
499 WRITE (KP, 8 1 5) 

GO TO 502 

....... FORMATS FOR INPUT AND OUTPUT STATEMENTS ...... 

10 FORMAT (4F 10.3) 

20 FORMAT (7F10. 5) 

30 FORMAT (2F 10. 3, 12, 8X, FI 0.2, 14) 

40 FORMAT (5F10.3) 

41 FORM AT (5A2,I5 , 6F 10. 5) 

42 FORMAT (6110) 

60 FORMAT (*1 9 ,40X, 9 PERFORMANCE OF A WIND-AXIS SIND TURBINE 9 ) 

61 FORMAT? 9 8 ,42X, 9 CALCULATED BY MODIFIED PROP PROGRAM 8 ) 

62 FORMAT (*0 9 ,8X,* BLADE SECTION - 9 , 5A2 , 5X , 9 NO. BLADES = %F2. 0 , 71, 9 R 

2ADIUS = 8 ,F5. 2, 8 METERS CONING ANGLE = *,F4. I, 9 DEGREES 8 ) 

63 FORMAT (* 0 9 , 8X, 9 METHOD OF CONTROL — *) 

64 FORMAT ( 8 *,3QX, 9 TIP SECTION PITCH, ANGLE OF PITCH = *,F5.1, S DEGREE 

4 S TIP SECTION LENGTH = 9 ,F4.1, 8 METERS 9 ) 

65 FORMAT? 9 * , 30X, 9 AI LERONS, SPAN = 9 ,F5.1, 9 METERS (OUTBOARD 8 , F5» 1 , 
5 9 PERCENT) DEFLECTION =*,F6.1, 9 DEGREES 9 ) 

66 FORMAT? 9 8 , 30 X, » SPOILERS, SPAN =’,F5.1, S METERS (OUTBOARD* , F5« 1 , 
6 8 PERCENT) DEFLECTION =*,F6»1, 9 DEGREES 9 ) 

67 PORMAT( 8 Q 9 ,35X, 9 AILERON CHORD IS 20% OF LOCAL BLADE CHORD 9 ) 

68 FORMAT ?* 0 9 ,30X, 5 SPOILER CHORD IS 10% OF LOCAL CHORD - - HINGED AT 
*70% CHORD*) 

69 FORMAT { 9 9 , 3X , 3 < 9 1 » ) I — • #5( « | -• ) , 

g 2 { 9 | —————— ® j f • | 3) 

70 FORMAT? 9 * # 3X l 3 ( 9 I — — — *) , 9 | — — -j— — *#3(* I -r— — « ) , 

71 FORMAT?* 9 ,3X, 9 | BLADE 1 RADIUS | LOCAL \ PITCH | ANGLE OFJ RES 
1 ULTANTf COEF. | COEF. 1 HINGE j TORQUE f A f ft** f*) 
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72 FORMAT (» '^ax,'] STATION 1 I CHORD 5 ANGLE ] ATTACK | VEL 

2 OCITY S OF 1 OF J MOMENT | INCREMENT) | S'? 

73 FORMAT { * »,3X,'| (PERCENT) | (METERS) j (METERS) | (DEG.) | (DEG.) | (H 

3/SEC) | LIFT 1 DRAG [ COEF. J (N-M) \ J ]*) 

74 FORMAT ( * ®,3X,'| FREESTREAM | ANGULAR ] TIP-SPEED] PITCH J TORQ 

40E 1 POWER J POWER J THRUST | TORQUE ) HINGE MOM ] A |®) 

75 FORMAT (® *,3X,'l WIND SPEED ] SPEED 1 RATIO | ANGLE S*#9X y 

5' | *,9X, • ) COEF. j',10X,«| COEF. J PER BLADE] (AVE) |*) 

76 FORMAT ( * »,3X,'] (MPH) ] (M/S) | (RPM) ] X | (DEG) | (N- 

6 M) | (KW) ] CP ] (NEWTON) ] CQ ] (N-M) J ]') 

77 FORMAT(»0' ,//,6X, 'INTEGRATED RESULTS --»,//) 

78 FORMAT (® 0* , 40X, ' NOTE; "A” IS ZERO ; INDUCED VELOCITIES ARE IGNORED* 

*//> 

79 FORMAT('0',40X, 'NOTE; -0.5.LE. A .LE.+0.5'//) 

80 FORMAT (' ' , F 1 0 . 1 ,F 10 „ 2, F 1 1 . 3, F8. 1 , F9 . 1 , FI 0. 1 , F9. 2 ,F 1 1 . 3 , F 1 0. 3 ,F 1 0. 
81 ,F9. 3,F9. 4) 

81 FORMAT ( * 0 * , F9. 1, F1 1. 2, F8. 1 , F10. 3 , F8. 1 ,F 1 2. 1 ,F9. 2,F 1 0. 3 ,F 12 . 1 ,F8. 3, 
1F10.3,F8.3) 

83 FORMAT(»0',40X,»NOTE:"A" IS NOT LIMITED'//) 

813 FORMAT (' 1 * ) 

815 FORM AT (' 0 * ,2X, * MODE' IN ERROR, PROGRAM STOPPED*) 

502 CONTINUF 
C 

STOP 

END 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
SUBROUTINE SEARCH (RL,C ,THET, CA,CS , D A, DS) 

SEARCH ~ DETERMINES THE CHORD AND THE TWIST ANGLE AT 

A GIVEN RADIUS ALONG THE SPAN. IT UTILIZES A LINEAR 
INTERPOLATION TECHNIQUE. 

COMMON /HOPE/ 1, DR,BB, B, V, X,THETP, AM OD, H, SI, GO, OMEGA, RHO, VIS, HL, 

1 PI, RX, W, NPROF , APF,RADEG,COSSI, CS2, NF,RR (50) ,CI(50) ,THETI(50) , 
2T1,T2,T3,T4,T5,T6,T7,T8,XETA,HH,B0,CAI (50) , CSI (50) , DAI (50) ,DSI(5Q) 
3 ,MCON,HMF,KFLAG 
RRV=RL 

IF (RRV. EQ. RR ( 1) ) GO TO 50 

LOCATE FIRST R IN ARRAY LESSTHAN RLOCAL, AND BRANCH TO 
INTERPOLATION STEPS 
DO 20 1=2, NF 

IF (RRV. GE. RR (I) ) GO TO 10 
C IF RLOCAL IS LESS THAN LAST VALUE IN TABLE, SET R=LAST VALUE 

IF (I. EQ. NF) GO TO 30 
20 CONTINUE 

10 J=I+1 

C COMPUTE INTERPOLATED VALUE 

PER= (RRV-RR (J-1) )/<RR{J-2) -RR (J-1) ) 

C=PER*(CI (J-2)-CI (J-1) ) +CI (J-1) 

THET=PER* (THETI (J-2) -THETI (J- 1) ) +THETI (J-1) 

CA=PER* (CAI (J— 2) — CAI (J- 1) ) +CAI (J-1) 

CS=PER* (CSI (J-2) -CSI (J-1) ) +CSI (J- 1) 

DA=PER* (DAI (J-2) -DAI (J-1)) +DAI (J-1) 

DS=PER* (DSI (J-2) -DSI (J-1)) +DSI(J-1) 
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RETURN 

C END INTERPOLATION 

30 C=CI (NF) 

C VALUES REQUIRED ARE AT FINAL STATION 

THET=THETI (NF) 

CA=CAI (NF) 

CS=CSI{NF) 

DA=DAI (NF) 

DS=DSI (NF) 

RETURN 

C VALUES REQUIRED ARE AT FIRST STATION 

50 C=CI(1) 

THET=THETI (1) 

CA=CAI ( 1 ) 

CS=CSI(1) 

DA=DAI (1) 

DS=DSI (1) 

RETURN 

END 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

SUBROUTINE CALC (RL,C,THET, QF ,TF , A # AP, F, CLF, CAT, CL, CD, ALFA, CH, CA, CS 

2, DA,DS) 

C 

C ... CALC - DETERMINES THE AXIAL AND ANGULAR INTERFERENCE 

C FACTORS AT A GIVEN RADIUS AND DETERMINES FUNCTIONS DEPENDENT 

C UPON THESE PARAMETERS. 

C 

COMMON /HOPE/ R, DR, HB, B,V, X,TBETP, AM OD, H, SI, GO, OMEGA, RHO, VIS, HL, 

1 PI,RX,W,NPROF,APF,RADEG,COSSI,CS2,NF,RR(50),CI(50) ,THETI (50) , 
2T1,T2,T3,T4,T5,T6,T7,T8, XETA, HH , BO ,C AI (50) ,CSI (50) ,DAI (50) ,DSI(50) 

3 , MCON , HMF , KFLAG 

C SET INITIAL VALUES 

DATA ACK/0.5/ 

XL=RL*OMEG A/V 
RH=HB 

SIG 8= B*C/PI/RL/8 „ 

THETA =THET+THETP 

IF LAST VALUE OF 'A* HAS GREATER THAN 0.5, SET INITIAL »A» EQUAL 
TO 0. TO START ITERATION. 

IF LAST VALUE OF *A» HAS EQUAL TO OR LESS THAN 0.5, USE LAST »A» 

TO START ITERATION 

BEGIN ITERATION FOR » A* AND *AP'. ITERATION WILL CEASE IF NOT 
CONVERGED AT 50TH STEP. 

IF (ABS ( A) .GT. ACK) A=0. 0 
DO 10 J=1 , 50 
C SAVE OLD 'A* AND A-PRIME (' AP ' ) 

BETA=A 

DELTA=AP 

C PHI IS RELATIVE HIND DIRECTION, INCLUDING INDUCED EFFECTS. 

PHI= ATAN ( (1.-A) *COSSI/ ( ( 1. +AP) *XL) ) 

IF (ABS (PHI) • LT. 0.0001) PHI= SIGN (0.000 1 , PHI) 

SINPHI=SIN (PHI) 

COSPHI=COS (PHI) 

C ALPHA IS ANGLE OF ATTACK 
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C THETA IS LOCAL BLADE PITCH ANGLE 

ALPHA=PHI- THETA 
C 

C ...... CALCULATION OF SECTIONAL LIFT AND DRAG COEFFICIENTS 

C 

C CALCULATE RESULTANT VELOCITY (' V 8 ) 

550 H=V*SQST ( C (1--A) *COSSX) **2+ < (1. *AP) *XL) **2) 

CALL NACAXX (R L, RX,SI , ALPHA, CL, CD, W) 

GO TO (22,2 1 #2 1) ,MCON 

C CALCULATE CONTROL SURFACE INCREMENTAL CL, INCREMENTAL CD, AND CH 

21 CALL INCHES (ALPHA,CL,CD ,CH, CA,DA, CS, DS,C) 

22 ALFA=ALPHA*57. 29578 
IF ( MCON. EQ. 1} CH=0, 

800 CONTINUE 

IFCGO.LT. 3.) GO TO 666 

CL=CLF*CL 

F=1 . 

GO TO 667 
C 

C ..... CALCULATION OF TIP AND HUB LOSSES ...... 

C 

666 IF<CAT. EQ. 1.) F=0.0 
IFCCAT.EQ. 1.) GO TO 667 
XXL= ABS (COSPHI/SINPHI) 

XXLO=XXL*R/RL 

CALL TIPLOS (XXL,XXLO,F, B,GO,HL,PI,R, RL,PHI,RH) 

C 

667 CX=CL*SINPHI-CD*COSPHI 
CY= CL*COSPHI+CD*SINPHI 

C * A MOD* DETERMINES MADEL FOR «A*. 

IF (A MOD. EQ. 0. ) GO TO 575 
¥BR=SIG8*CY*CS2/ (SINPHI**2) 

¥AR=SXG8*CX/F/SINPHI/COSPHI 
CAN=F*F+4. *VBR*F* (1.-F) 

I F (CAN. LT, 0.0) CAN=0. 0 

C S A» IS INDUCED AXIAL VELOCITY FACTOR, *A» IS INDUCED RADIAL 

C VELOCITY FACTOR. 

A=C2.*VBR + F-SQRT (CAN) ) / (2. * (VBR+F*F) ) 

C I? 1 APF 9 IS 1.0, PROGRAM CONVERGES ON 8 A*, NOT 8 AP 8 « 

C IF 8 APF * IS NOT 1.0, PROGRAM CONVERGES ON 8 AP 8 , NOT «A». 

IF (APF. EQ. 1.) GO TO 580 
AP=VA R/ ( 1 » -VAR) 

GO TO 580 

575 VBR=SIG8*CY*CS2 
VAR=SIG8*CX 
IF I KFLAG. EQ. 1) A=0.0 
IF (KFLAG. EQ. 1 ) GO TO 579 
A=VBR/ (F*SINPHI**2+VBR) 

579 IF (APF. EQ. 1.) GO TO 580 

AP=VAB/(F*SINPHI*COSPHI-VAR) 

580 CONTINUE 

IF (KFLAG. EQ. 3) GO TO 581 
IF{A.GT.0.5) A=0. 5 
IFfA.LT. (-0.5) ) A=-0.5 
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C ....... DAMPENING OF AXIAL AND ANGULAR INTERFERENCE FACTOR 

C ITERATIONS* 

581 CONTINUE 
40 A= (A+BETA) *«5 

AP= (AP+DELTA) *.5 
30 CONTINUE 

C 

C ....... TEST FOR CONVERGENCE ...... 

C 

I F ( APF. EQo 1.) GO TO 70 
IF (ABS(AP-DELTA) .LE..001) GO TO 50 
GO TO 10 

70 IF(ABS{A-BETA).LE..001) GO TO 50 

C WHEN SPECIAL DIAGNOSTIC PRINT-OUT IS NEEDED, CARDS ARE INSERTED HE 

C 

10 CONTINUE 

KP=6 

WRITE (KP,79) 

79 FORMAT* *0* ,2X,' VALUE OF SLOWDOWN FACTOR, A , IS NOT CONVERGED*/) 

50 W=V*SQRT (( (1.-A) *COSSI) **2+ ( ( 1. +AP) * XL) **2) 

CT1= {0. 5*RHO*B*C) * (W*W) 

QF=CT1*RL*CX 
TF=CT 1*CY*COSSI 
HMF=CH*CT1*CA**2/C 
C 

RETURN 

END 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
SUBROUTINE TIPLOS <U,UO,F,Q» GO, HL, PI, R,RL,PHI,RH) 

C 

C ...... TIPLOS - DETERMINES THE TIP AND HUB LOSSES 

C BASED UPON GOLDSTEIN'S THEORY, OR PRANDTL'S THEORY, 

C OR FOR THE CASE OF NO LOSSES. 

C 

SUM2=0.0 
SUM=0. 0 
AK=1. 

AMM=1. 

AM=0. 0 

IF (Q.GT. 2.0) GO TO 966 
IF(GO.EQ.O.O) GO TO 200 
IF (GO.EQ. 1.0) GO TO 100 
IF {GO. EQ. 2. 0) GO TO 444 
966 IF (GO. EQ. 2.0) GO TO 444 
200 CONTINUE 

F= (2. /PI) ♦ARCOS (EXP (- (Q* (R-RL) ) / (2. *RL*SQRT (SIN (PHI) **2 +. 000 1) ) ) ) 
GO TO 105 
444 F=1 .0 

GO TO 105 

100 IF ( (ABS (SIN (PHI) ) ) . LT. . 000 1 ) GO TO 200 
C 

C .... GOLDSTEINS METHOD....... 

C 

DO 10 M=1 , 3 
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V=(2.*AM+1.) 

ZO=(IO*V 

V2=V*¥ 

Z=0*V 

Z2=Z*Z 

CALL BESSEL (Z,V,AI) 

CALL BESSEL (ZO*V,AIO) 

IF (Z. GE. 3. 5) GO TO 300 
A=2.*2. 

B=4.*4. 

C=6.*6. 

D=8. *8, 

T1¥Z=Z2/(A-V2) + (Z2*Z2)/((A-?2) *(B-V2) ) + (Z2**3)/( (A-V2) * (B-V2)* 

1 (C-V2)) + {Z2**4)/((A-V2)*(B-V2)*(C-V2)*(D-V2) ) 
CT1VZ=(Y*PI*AI)/(2.*SIN (. 5*V*PI) ) -T1YZ 
GO TO 400 

300 T0=(0*0)/(1.+0*0) 

T2=4. *0*0* (1 .-0*0) / ( (1. +0*0) **4) 

T4=16. *0*0* (1.-14. *0*U + 21. *0**4 - 4. *0**6)/ ( (1. + 0*0) **7) 

T6=64. *0*0* (1.-75. *0*0+603. *0**4-1 06 5. *0**6 +4 60. *0**8- 36. *0**10) 
1/ ( ( 1 - +-0*0) **10) 

CT1VZ=T0+T2/V2+T4/(V2**2)+T6/(V2**3) 

400 FV0=(0*0)/(1.+0*0) - CT1VZ 

SOM=SOM+FVO/V2 
IF (AM. HE. 0.0) GO TO 1 
E=-0. 098/ (00**. 668) 

1 IF(AM.NE. 1.0) GO TO 2 
E=0. 031/(00**1. 285) 

2 IF (AM. GT. 1.0) E=0.0 

SOM2=SOH2+ { <00*00*AHM)/(1. + U0*00) - E)*(AI/AI0) 

AM=AM+1 « 

AK= ( (2. *AM~ 1.)*AK)/(2.*AM) 

10 AMM=AK/(2.*AM+1.) 

G= (0*0) / (I. *0*0) - (8./ (PI*PI) ) *SOM 
CIRC=G— (2. /PI) *S0M2 
F= ( (1. +0*0) /(0*0) ) *CIRC 
C 

C HOBLOSS CALCOLATIONS 

C 

105 IF(HL.EQ. 1.0) GO TO 500 
FI= 1.0 
GO TO 900 

50 0 FI= (2./PI) *A8C0S (EXP (- (Q* (RL-RH) ) / (2 . *RH*SQRT (SIN (PHI) **2+ 

1 . 0001 )))) 

900 F=F*FI 

RETORN 
END 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccCcccccc 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

SOBROOTINE BESSEL (Z, V, AI) 

c 

C BESSEL CALCOLATES BESSEL FONCTIONS FOR THE GOLDSTEIN 

C TIP LOSS MODEL 

C 
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S=0.0 
AK=0 - 0 
C=1. 

DO 30 K= 1 f 10 
B={. 25*Z*Z) **AK 
D=V+AK 
P=1. 

5 TK=D~ 1 , 

IF (TK.LE. 0.0) GO TO 40 

p=D*TK*P 

D=D-2. 

GO TO 5 
40 E=P 

S=B/ (C*E) + S 
AK=AK+1. 

C=AK*C 

30 CONTINOE 

AI= ( (. 5*Z) **V) *S 
RETORN 
END 

SOBROOTINE NACAXX (RL, RX ,SI , ALPHA, CL, CD, W) 

THIS SOBROOTINE CALCOLATES LIFT AND DRAG COEFFICIENTS FOR THE NACA 
23024 HALF ROOGH AIRFOIL 

AS A FONCTION OF THE ANGLE OF ATTACK ALPHA (IN RADIANS). 

Y = ALPHA* 180./3. 14 15927 
AALPHA = ABS (ALPHA) 

X = ABS (Y) 

CONTINOE 

IF (8 . .LT. X .AND. X ,LT.13.)CD = X/248. 8-. 0 1675 

IF (-9. .LE. Y .AND. Y -LE. +8.) CD = -2. 09050E-6*Y**3 4 - 7 . 19578E-5*Y 
***2 +2. 26766E”4*X +9.99331E-3 

IF f 1 3. .LE. X .AND. X .LT. 160.) CD = 2. - 1. 089* (1. 57- AALPHA) **2 
IF ( 1 60. .LE. X .AND. X .LE. 180. ) CD = 0.04 + 2.8* (AALPHA-3. 142) **2 


CONTINOE 

IF ( 1 0. -LE- X .AND. X . LE. 16.)CL * -1. 25E-3*X**3 +4. 55953E-2*X**2 
* -5. 47084E- 1*X +3.04822 
C 

IF (~12. .LE. Y .AND. Y .LT. 10.) CL = 2 . 646 1 9E- 10* Y**8 +2.39782E-8* 
* Y**7 +1 » 03975E-7 *Y**6 

1~5. 54091E— 6*Y**5 -2.20717E-5*Y**4 +3 . 1 1 34 1E~4*Y**3 +4. 591 1 0E-4*Y** 
*2 +8.26369E-2*v +1.00579E-1 
C 

IF (16. .LT. X .AND. X . LE. 90.)CL = 1.1- 1 . 7 8* (AALPHA- . 7853) **2 
C 

IF(90. .LT. X .AND. X . LE. 160. ) CL = -1.1+ 1. 78* (AALPHA-2. 356) **2 
C 

IF (160. .LT. X .AND. X .LE. 172.5)CL = -.763 
C 
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IF ( 172. 5 . LT. X .AND. X „ LE„ 180. } CL = .10173*1 - 18.3114 
C 

IF (Y .LT. ~ 12. } CL = -CL 
C 

RETURN 

END 

CCCCCCCCCCCCCCCCC SUBROUTINE ISCREM CCCCCCCCCCCCCCCCCCCC 
C NRITTEN BY ». H. NENTZ AND CYRUS OSTAHARI, HSU, JULY 1979... 

C ADJUSTS AIRFOIL COEFFICIENTS FOR DEFLECTION OF SPOILERS OR AILER 

SUBROUTINE INCREM{ALPHA,CL,CD,CH,CA, DA,CS,DS,C> 

PI=3. 1.41593 
Y= ALPHA* 1 80. /PI 

C CHO DEPENDS ON AIRFOIL CAMBER NEAR TRAILING EDGE. 

C THIS VALUE OF CHO IS FOR NACA 230 MEAN LINE AILERONS ONLY! 

CHO=-. 0035 
CH=0. 

DAR=DA*PI/180. 

DSR=DS*PI/180. 

IF (CA . EQ. 0. ) GO TO 5 

C BEGIN AILERON INCREMENT CALCULATIONS 

C CHECK FOR POST STALL CONDITION 

IF (17 . . LT. Y.OR. Y.LT.- 17. ) GO TO 25 
IF (0. • LE. DA. AND. DA. LE.60.)GO TO 10 
IF(-60. . LE. DA. AND. DA.LT.0.) GO TO 15 
10 DCL=CA/C* (11. 5*SIN (DAR) -5. *SIN {D AS) * SIN (DAR) ) 

DCD=1.05*CA/C*SIN (DAR)*SIN (DAR) 

CH=CA/C* (CHO— . 0542* Y- 2. *SI N (DAR) ) 

GO TO 20 

15 DCL=CA/C* ( 1 1. 5*SIN (DAR) *3. 35*SIN (DAR) *SIN (DAR) ) 

DCD=. 65*CA/C*SIN (DAR) *SIN (DAR) 

CH=CA/C* (CHO-.0542*Y-2.*SIN (DAR) ) 

GO TO 20 

5 IF (CS. EQ. 0. ) GO TO 30 

C BEGIN SPOILER INCREMENT CALCULATIONS 

IF(17.LT. Y.OR. Y.LT. -17. )GO TO 30 
C SPOILER IS ASSUMED INEFFECTIVE FOR POST-STALL CASES 

DCL=-9.9*CS/C*SIN (DSR) 

DCD=. 93*CS/C*SIN (DSR) *SIN (DSR) 

IF(0..LT. Y.AND.Y.LE. 17. ) CH=- ( 6. 7-. 17 5*Y) *CS/C*SIN (DSR) *SIN (DSR) 
IF(-17..LE.Y. AND. Y.LE.O.) CH=- (6. 7-, 075*1) *CS/C*SIN (DSR) *SIN (DSR) 

20 CL=CL»DCL 

CD=CD+DCD 
GO TO 30 
25 CONTINUE 

C CALCULATE POST-STALL CONDITIONS FOR AILERON 

C CONVERT CL AND CD TO CN AND CC 

C 

C NORMAL FORCE COEFFICIENT: 

C 

CN=CL*COS (ALPHA) +CD*SIN (ALPHA) 

C 

C CHORDHISE FORCE COEFFICIENT: 

C 

CC=CD*COS (ALPHA) -CL*SIN (ALPHA) 
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CALCULATE CONTRIBUTIONS TO FORWARD ELEHENT 

CM { 1 „-~CA/C) 

CC1=€C* (1 »-€&/€) 

CALCULATE AILERON CONTRIBUTIONS 

CN2=2»27*2»*PI*SIN(ALPHA+DAR) /{4,4-PI *ABS(SIN (ALPHA^DAR) ) ) 
DCNA=CN2*CA/C*C0S (DAR) 

DCCA=CN2*CA/C*SIH (DAR) 

CH=~CN2/2« 

CN=CMl-s-DCNA 
CC=CC1*DCCA 

CONVERT TO LIFT AND DRAG COEFFICIENTS: 

CL=CN*CQS (ALPHA) -CC*SIN (ALPHA) 

CD=CC*COS (ALPHA) +CN*SIN (ALPHA) 

30 CONTINUE 
RETURN 
END 
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Cb) Incremental Drag 


Figure A- 2 - Continued, 
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O 12° 

bv 16 ° 


( c) Hinge Moment 
Figure A- 2- Concluded. 
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Figure A- 3 - 10% Spoiler Incremental Effects. 








APPENDIX B 


SYMBOLS 


a axial induced velocity factor 

a' rotational induced velocity factor 

A hydraulic piston area, sq in 

A r aileron rib area, sq in 

A s aileron skin area, sq in 

B number of blades 

c airfoil or blade local section chord, m (ft) 
c a aileron chord, m (ft) 

c^ hinge moment coefficient 

c^q hinge moment coefficient at a = 0° and 6=0° (a function of camber) 

Cfo rate of change of hinge moment with alpha, obtained from the data of 

a reference 5 

c h „ rate of change of hinge moment with aileron deflection, from 
reference 5 

ca airfoil lift coefficient 

c n normal force coefficient 

c s spoiler chord, m (ft) 

c s airfoil chordwise force coefficient 

c x coefficient of force component in plane of rotation 

Cy coefficient of force component in axial direction 

Cp rotor power coefficient, powerAp V 3 tt R 2 

D hydraulic piston diameter, in 

F actuating cylinder force, lb 

F c total centrifugal force on aileron, lb 
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g 


acceleration of gravity, 32.2 ft/sec 2 
H airload hinge moment, in-lb 

H c hinge moment supplied by centrifugal force, in-lb 
K correction factor, Kirchoff flow 

L length of hydraulic piston travel, in 
m aileron mass, slugs 

n load factor = 

y 

N aerodynamic normal force, N (lb) 

p hydraulic pressure, psi 

q jp (wx) 2 , dynamic pressure due to blade motion, N/m 2 (lb/ft 2 ) 
r local blade radius, m (ft) 

R maximum blade radius, m (ft) 

t skin thickness, in 

T aileron deflection torque, in-lb 

V wind velocity, m/sec (mph) 

w weight, lb; also specific weight of aluminum, lb/in 3 

w r weight of ribs, lb 

w s weight of skin, 1b 

x spanwise distance from center of rotor, m (ft) 

X tip speed ratio, u>R/V 

a angle of attack, deg 

<5 deflection angle, deg 

<j> relative wind angle, deg 

e blade pitch angle, deg 

(o rotor angular velocity, rad/sec (rpm) 
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